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ABSTRACT 

In a recent study (Martizzi et al. 2012), we used cosmological simulations to show that 
active galactic nuclei (AGN) feedback on the gas distribution in clusters of galaxies can 
be important in determining the spatial distribution of stars and dark matter in the 
central regions of these systems. The hierarchical assembly of dark matter, baryons and 
black holes obscures the physical mechanism behind the restructuring process. Here 
we use idealized simulations to follow the response of a massive dark matter halo as 
we feed the central black hole with a controlled supply of cold gas. This removes most 
of the complexity taking place in the cosmological simulations that may have biased 
our previous study. We confirm our previous results: gas heated and expelled from the 
central regions of the halo by AGN feedback can return after cooling; repeated cycles 
generate gravitational potential fluctuations responsible for irreversible modifications 
of the dark matter mass profile. The main result is the expulsion of large amounts 
of baryons and dark matter from the central regions of the halo. A strong prediction 
of this scenario is similar sized cores in both the central stellar and dark matter 
components of massive galaxies. 

Key words: black hole physics - cosmology: theory - cosmology: large-scale structure 
of Universe - galaxies: formation - galaxies: clusters: general - methods: numerical 



1 INTRODUCTION 

In the concordance cosmological scenario, the large scale 
structure of the universe is globally determined by the dom- 
inant mass component, dark matter. The situation is differ- 
ent within highly non-linear regions where the baryons dom- 
inate. Gas cooling, star formation, feedback processes and 
dynamical interactions between different mass components 
can modify the spatial distributions of all components. For 
example, dark matter halos are known to globally respond 
to the condensation of baryonic matter in their central re- 
gions via dissipative processes. This increases the central 
mass density through adiabatic contraction, and causes the 
halo to become rounder. More recent studies of baryonic 
processes have f ocussed on the effects of supernovae feed- 
back in galaxies (iGovernato et al]|2010l ; iMaccio ' et al.ll201ll : 
iTevssier et al.ll2012l h however there are fewer studies of the 
effects of AGN feedback on massive galaxies. 

Clusters of galaxies provides a wealth of opportunities 
to test the validity of the current paradigm for cosmic struc- 
ture formation. Analysing the details of the spatial distribu- 
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tion of the different components of a cluster, as predicted by 
theory, and comparing them to those measured (or inferred) 
from observations allows us to test our current understand- 
ing of high mass galaxy formation. 

In this paper we focus on the role of active galac- 
tic nucleus (AGN) feedback in shaping the mass distribu- 
tion of clusters of galaxies. The main motivation for this 
work comes from recent observational c onstraints of cluste r 
mass profile s using gravitational len s ing: | Sand et al.l (2004), 
ISand et ail ll200ah [Newman et~al] (120091'), [Newman et al l 
j|201ll ), iRichtler et all (|201lh . iNewman etafl (|2012[ ). These 
report mass density profiles whose central slopes, measured 
within the inner ~ 5 kpc, are much shallower than the 
universally adopted NFW profile. On a slightly smaller 
scale, within 1-2 kpc from the centre, very shallow slopes 
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lGrahamll201lf ). The combination of these two observational 
facts suggests a possible connection between the mechanisms 
that produce shallow density slopes in the dark matter dis- 
tribution and in the stellar distribution. For example, both 
the dark matter and stellar components can respond simi- 
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larly to dynamical perturbations since they behave like col- 

lisionless fluids. 

In a recent study (jMartizzi et ahlfeoid ). we used high 
resolution cosmological hydrodynamical simulations includ- 
ing AGN feedback to show that it is possible to produce 
clusters with similar properties as indicated by these ob- 
servations. We argued that at least three processes could 
contribute to creating shallow density profiles: (i) super- 
massive black holes (SMBHs) can transfer part of their 
orbital energy to collisionless matter via dynamical fric- 
tion during mergers; part of the coll isionless matt e r is e x- 
pelled from the central regions (i.e. iGoerdt et ail l |201Ch ); 
(ii) AGN driven gas outflows modify the gravitational 
potential close to SMBHs; ejection of collisionless mat- 
ter from the central region of the cluster is produced if 
the potential fluctuations happen on a time-scale smaller 
than, or resonant with, the dynamical time. To be effec- 
tive this me chanism requi r es many subsequent p o tentia l 
fluctuations dNavarro et al.l il 19961); iGnedin fc Zhaol d2002h : 



iRead fc Gilmorel l |2005h ; |Pontzen fc Governatol l |2012h ); (hi) 
slow loss of gas mass from the central regions of a cluster 
can produce an expansion of the inner mass distribution; 
this slow loss of gas can be obtained as an effect of AGN 
heating. 

Our previous results indicated that process (ii) was 
the dominant contribution to the shallowing of our cen- 
tral dark matter and st ellar cusps. Very similar re- 
sults were also found by iDuffv et all (|20ld ) using the 
OWLS suite of cosmological simulations. Related mech- 
anisms have been found t o be active in lower mass 
halos (Governato et al.l |2010|; iPontzen ~Ei Governatol |2012t 
iMaccio' et all l201ll ; ITevssier et alj|2012h . where the role of 
AGN feedback is substituted by that of supernovae explo- 
sions. Unfortunately, the complexity of hierarchical struc- 
ture formation makes it difficult to disentangle the efficiency 
of this mechanism. To obtain a clearer understanding we will 
use simpler, idealised simulations, that remove these ambi- 
guities in order to provide a better insight into the core 
creation process. 

In a recent paper, iRagone-Figueroa et al.l (|2012l ) study 
the role of baryonic mass outflows in galaxy clusters, us- 
ing collisionless matter simulations where gas is modeled 
as a time varying external potential. In this paper we 
choose a similar approach, but we choose to perform full N- 
body+hydrodynamical simulations of a isolated dark matter 
halo, including gas cooling and AGN feedback from a single 
central supermassive black hole. AGN heating can produce 
the impulsive expulsion of a large quantity of gas from the 
central regions of a cluster. The result of this is a fast varying 
central gravitational potential without the dynamical effects 
of hierarchical mergers of galaxies and their black holes. 

The paper is organized as follows: the first section is 
dedicated to the numerical methods we adopted for our sim- 
ulations; the second section gives our the main results; in the 
last section we discuss and provide a short summary. 



2 THE SIMULATIONS 

We use the AMR code RAMSES (|Tevssierl [20021 ') to per- 
form simulations that follow the evolution of the mass dis- 
tribution of an isolated dark matter halo in presence of a 



gaseous component. The simulations differ in the way ini- 
tial conditions for the gas are set, as explained in subsection 
12.11 The gas dynamics i s modeled using a second-order un- 
split Godunov schem e (|Tevssierl 120021 : ITevssier et al.l 120061 : 
iFromang et al.|[2006h based o n the HLLC Rie mann solver 
and the MinMod slope limiter (|Toro et al.ll 1994 ) . We assume 
a perfect gas equation of state (EOS) with 7 = 5/3. We sam- 
ple the dark matter distribution with 10 6 particles of mass 
m-dark = 1-78 x 10 7 Mg> . In all simulations gas represents 
17% of the total mass in the halo; the gas mass resolution 
element is m gas = 3.79 x 10 6 M©. The minimum cell size is 
Axmin = 0.18 kpc. The AMR grid was dynamically refined 
using a quasi-Lagrangian strategy: when the dark matter or 
baryons mass in a cell reaches 8 times the mass resolution, 

it is split into 8 children cells. 

Gas is allowed to cool: we use the lSutherland fc Dopital 

Il993l cooling function to account for H, He and metal cool- 
ing. We include AGN feedback, using the s ame implemen- 
tation we adopted in Tevssier et al.l (|201 lh . iMartizzi et al.l 
(120121) and IMartizzi et"aL 

I (120121 1. a modified version of the 
iBooth fc Schavd i2009l) model. SMBHs are modeled as sink 
particles, following iKrumholz et al.l (j200J). The gas accre- 
tion onto SMB Hs is computed using a modified Bon di-Hoyle 
formula jBooth fc Schavell2009l : [Dubois et al.ll2012l ). A frac- 
tion of the accreted mass is converted into thermal energy 
that is directly injected into the gas surrounding the black 
hole. Despite its simplicity, when properly tuned, the model 
reproduces basic features of AGN feedback (as it is thought 
to work), in particular the existence of two feedback modes: 
the impulsive 'quasar mode' during cold gas accretion when 
the accretion rate is high, and the quiescent 'radio mode' 
during hot gas accretion when the accretion rate is low. In 
this picture, the effect of 'radio mode' feedback is to prevent 
hot gas from cooling, while the effect of 'quasar mode' feed- 
back is to eject gas from galaxies, violently suppressing star 
formation. 



2.1 Initial conditions 

The initial conditions for dark matter are sca led-up ver- 
sions of those adopted by ITevssier et al.l (|2012T ) . The dark 
matter distribution follows the NFW density profile. We 
chose a concentration parameter c = 10 and a circular 
velocity V200 = 350 km/s, corresponding to a virial mass 
M200 = 1.42 x 10 13 M Q and a virial radius i? 2 oo = 500 kpc. 
We truncate the halo at 1.13 Mpc, so that the total enclosed 
mass is Mhaio = 1.78 x 10 13 Mg. Progenitors of halos of mass 
~ 10 14 Mq identified at redshift z = have a typical mass 
of a few 10 13 M at redshift z > 2. This means that our 
simulations are particularly suited for comparing with the 
early evolution of clusters of galaxies. We choose the initial 
velocity field of the particles in a way that equilibrium can 
be maintained and such that the system has spin parameter 
A = (no rotation). This choice has been made to avoid 
dealing with the details of gas acc retion onto SMBHs in a 
clumpy disk (|Bournaud et al.ll2012l ) and it allows us to study 
test cases in which a high gas accretion efficiency is reached. 

In all simulations a gaseous halo is embedded in the 
dark matter halo. The gas density p, pressure P and tern- 
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perature T follow specific profiles: 



p(x) = pa 



ln(l + x) 



P{x) 
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T(x) = T 



ln(f + x) 
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where r s is the NFW scale radius, p s is the dark matter 
density at r s , and x = r/r s . We set the central temperature 
To according to 



k B T 2 T 

= 4nGp a r s — 

pm p 
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(4) 



Note that our choice implies that P oc p . The effective 
polytropic index F is set according to 



f + 



(1 



)ln(l 



(5) 



(f + 3x eq ) ln(l + x eq ) 

where x eq = \/5c. See lKomatsu fc Seliakl (120011 ) for a sim- 
ilar, but slightly more complex model. With our choice of 
parameters, the gas density profile follows the slope of the 
NFW profile for x > x eq and tends to a constant value 
po for x <C x eq . The central density po is set requiring the 
gaseous halo mass fraction to be /g aa = 0.15. The metallic- 
ity of the gaseous halo is set to Z — 10~ 2 Zq. The stability 
of the initial conditions was confirmed by running adiabatic 
simulations (no gas cooling). 

Within the halo we place ten spherical gas clumps of 
equal mass, AI C \ = 2.1 X 10 10 Mq and radius R c \ — 10 kpc. 
We set the temperature of the gas in the clumps to be a few 
10 4 K; its metallicity is set to Z = Zq. Each clump is given 
an initial velocity of magnitude V c i = V200 directed towards 
the halo center. The aim of these clumps is to roughly mimic 
the cold gas reservoirs provided by infalling galaxies into a 
cluster. We have chosen radial orbits to be able to efficiently 
provide gas to the central regions of the halo; such a choice 
allows us to study the maximum effect gas dynamics can 
have on the dark matter distribution. 

Our simulations differ in the initial geometrical distri- 
bution of the clumps. In one simulation, the distance of the 
clouds from the centre is in the range 200 kpc < r < 1 Mpc; 
in this case the clouds far from the centre and reach it slowly, 
we labeled this simulation as SLOW. In the second simula- 
tion the distance of the clouds from the center is in the range 
50 kpc < r < 250 kpc; since the clouds are close to the cen- 
tre and reach it faster, we labeled this simulation as FAST. 
We also consider a third run in which the clumps distribu- 
tion is the same as the SLOW case, but no AGN feedback 
is considered; this run has been labeled as AGNOFF. 

Finally, as additional test cases we run two simulations 
which only include the dark matter and smooth gaseous ha- 
los. The first one includes AGN feedback and is labeled as 
SMOOTH AGNON, the second one does not include AGN 
feedback and is labeled as SMOOTH NOAGN. In these two 
runs the halo metallicity is set to a higher value Z — Q.2Zq. 
These latter two simulations allow us to study the effects of 
a cooling flow, rather than clumpy accretion. 

In the FAST, SLOW and SMOOTH AGNON simula- 
tions, we place a SMBH of mass Mbh = 10 9 Mq at the 
centre of the halo. This choice for the black hole mass 



roughly matches the value measured at redshift z ~ 2 for 
th e central SMBH in t he cosmological simulation analysed 
in lMartizzi et all (|2012T l and agrees with the expe cted mass 
from the Mbh — o relation (|Tremaine et al.ll2002r i. 



3 RESULTS 

We run our simulations for 3.5 Gyr. For comparison, the 
dynamical time in the central region of the halo (r < 5 kpc), 
which depends on the density and radius, is of the order of 
10 s yr. As shown by the top left panel of Figure [T] 3.5 Gyr 
is sufficient time to see convergence in the effects of dark 
matter removal from the central regions of the halo. The 
black continuous line shows the initial mass profile in both 
simulations. After 1-2 Gyr a significant amount of mass has 
been ejected from the region r < 5 kpc in both the FAST 
(red line) and SLOW (green line) simulations. 

The convergence to a central inner core occurs in both 
simulations with AGN feedback. On the contrary, the AG- 
NOFF simulation (blue line) is completely different: as the 
clumps fall to the central region of the cluster and as the 
gaseous halo cools, the dark matter halo contracts as a re- 
sponse to gas condensation and as a result the central den- 
sity increases. 

Some differences can be observed when considering the 
evolution of the mass of the central SMBH as a function 
of time t in our models (top right panel of Figure [T]). As 
expected, in the FAST simulation the SMBH starts grow- 
ing much earlier. The final black hole masses are Mbh = 
2.28 x 10 9 M Q and Mbh = 2.86 x 10 9 M for the FAST 
and SLOW simulations respectively. In the FAST simulation 
the gas clumps feed the SMBH earlier, triggering a phase of 
rapid mass growth (t < 1.5 Gyr) followed by a phase of in- 
tense AGN feedback. This is then responsible for the subse- 
quent slowing down of mass accretion onto the SMBH. In the 
SLOW simulation, the gas clumps start feeding the SMBH 
later, so the early evolution of its mass is mild (t < 1.8 Gyr); 
at later times a large quantity of gas becomes available to ac- 
crete onto the SMBH, triggering a very rapid mass growth 
in the time interval 1.8 Gyr < t < 2.1 Gyr; this phase is 
also followed by intense AGN activity which slows down the 
SMBH mass growth. 

The bottom panels of Figure[T]show the evolution of the 
gas and dark mass enclosed within 5 kpc from the halo cen- 
tre as a function of time. This is useful for understanding 
the effect of AGN feeback in the central region. The pas- 
sages of the gas clumps close to the centre during the early 
evolution can be detected as peaks in the enclosed gas mass 
(t < 0.7 Gyr for FAST and t < 1.4 Gyr for SLOW). Before 
the central SMBH becomes massive enough to unbind the 
gas provided by the clumps, a large quantity of mass is accu- 
mulated in the central region. This produces a contraction 
of the total mass distribution that can be observed in this 
plot as a steep increase in the gas and dark mass within 5 
kpc. 

The striking effect we observe is that this contraction 
is completely erased as soon as the AGN activity becomes 
violent enough to unbind the gas. After the first passage of 
the clumps through the centre, the high metallicity material 
mixes with the low metallicity gas in the halo. By this time, 
a large reservoir of cooling gas is available. This gas is feed- 
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Figure 1. Top left: dark matter density profiles at different times and for different simulations. After 3.5 Gyr a large amount of mass 
has been removed from the central region in the FAST and SLOW simulations. The AGNOFF simulation shows signs of adiabatic 
contraction, with a large accumulation of mass in the central region. The solid lines are obtained when the peak of the dark matter 
distribution is chosen as a center, whilst the dashed lines are obtained when the SMBH is chosen as the center. Top right: SMBH mass 
evolution in the FAST and SLOW simulations. The mass evolves differently in the two runs, but the qualitative evolution through bursts 
of Eddington limited gas accretion is the same. Bottom panels: evolution of the dark (red line) and gaseous (green line) mass within 5 
kpc from the cluster center in the FAST (left) and SLOW (right) simulations. After 2 Gyr a large amount of dark matter is removed 
from the central region. The inset panels show a zoomed view of the evolution of the gas mass within 5 kpc, which correlates with the 
evolution of the accretion rate. As soon as gas is pushed away from the SMBH by AGN feedback the accretion rate drops. 



ing the SMBH but is also influenced by its feedback. AGN 
heating results in the ejection of gas from this region. The 
ejected gas eventually cools down and begins infalling again. 
The accretion-ejection cycle happens on a timescale that is 
always comparable to the dynamical time in the considered 
region. Fluctuations of the dark matter mass enclosed in this 



region follow the fluctuations in the gas mass. The result of 
subsequent accretion-ejection cycles is the irreversible ejec- 
tion of ~ 2 x 10 10 Mq from r < 5 kpc, which is roughly 60 
% of the dark mass initially enclosed in this region. 

There are interesting differences in the properties of the 
accretion-ejection cycles in the two simulations: they are 
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Figure 2. Top panel: distance of the SMBH from the centre of the 
dark matter distribution as a function of time. Bottom: gas mass 
within 5 kpc as a function of time. The SMBH is chosen as the 
centre. This plot shows that the decrease in the gas mass poorly 
correlates with the black hole motions. Therefore, the decrease 
in gas mass is not a consequence of the SMBH moving to lower 
density regions but rather due to AGN feedback. 



longer and more irregular in the FAST case, but they ap- 
pear to be quasi-periodic in the SLOW case. This difference 
is possibly due to the fact that for t > 2 Gyr the SMBH in 
the SLOW simulation can provide a more efficient feedback 
due to its larger mass. 

The inset panels in Figure Q] show the correlation be- 
tween the accretion rate onto the SMBH and the amount 
of gas in the centre. Whenever gas mass in the centre is 
removed, the accretion rate suddenly decreases. To check 
whether this is a real effect or just a product of the black 
hole oscillating within a perturbed potential or within a har- 
monic core, we compare the distance of the SMBH from the 
dark matter centre as a function of time to the evolution 
of the gas mass within 5 kpc during the same time interval 
(Figure [5]). This plot is made using the SLOW simulation 
and a similar result is found in the FAST case. There is no 
obvious correlation between these quantities, implying that 
the adopted centering technique is not biasing or results. 

Figure [3] shows a time sequence of images of the gas 
distribution over an interval of 100 million years. This spans 
the period in which the SMBH dramatically lowers the cen- 
tral gas density. The left panels show four projected gas 
density maps at four different times, before, during and af- 
ter an intense AGN activity burst. The SMBH is located 
at the centre of each map. The right panel of the same fig- 
ure shows the evolution of the gas mass within 5 kpc (green) 
from the SMBH and the re-scaled gas accretion rate (black). 
The times of the above four snapshots in the maps are la- 
beled as numbers in the right panel. We can clearly see that 
a large amount of gas is available for accretion in snapshot 



1 and this is observed in the right-hand plot as large values 
of the accretion rate and enclosed gas mass within 5 kpc. In 
snapshots 2 and 3, AGN feedback blows gas away from the 
central region and the accretion rate decreases. Finally, in 
snapshot 4 we see gas returning to the centre and beginning 
to accrete again onto the SMBH. 

In Figure 2] we plot a = dlog(p)/dlog(r), the central 
slope of the dark matter density profiles as a function of 
time. We fit a power law to the density profiles in the region 
0.4 kpc < r < 8 kpc, comparing the central slopes we obtain 
by choosing the dark matter peak as the centre (left panel) 
to those we get by choosing the SMBH at the centre (right 
panel). The dark matter peak location has been measured 
using a shrinking spheres algorithm. As it is well known, in- 
ferred density profiles are sensitive to the choice of centre. 
When the SMBH is chosen as centre, the slope is approx- 
imately a = —0.2, close to a constant density. A slope of 
a = —0.7 is found if we choose the dark matter peak as the 
centre. 

Both the FAST and SLOW simulations converge to the 
same final density profiles and their evolution over time is 
qualitatively similar. Initially, the clumps transfer their en- 
ergy to dark matter via dynamical friction; the result of this 
energy transfer is an initial shallowing of the slope of the 
dark matter profile. This p r ocess has been shown to be effi- 
cient by e.g. lEl-Zant et all \2004 ). As gas in the halo cools, 
this adiabatically contracts the dark matter and steepens the 
central profile, reversing the initial flattening. This period 
continues until the black hole has grown massive enough to 
halt the accretion process. The steepest slopes are obtained 
at t ~ 0.6 Gyr and t ~ 1.2 Gyr in the FAST and SLOW 
simulations respectively, as expected from the difference be- 
tween the clump arrival times in the two runs. 

As soon as AGN feedback becomes strong enough to 
unbind large quantities of gas, the accretion-ejection cy- 
cles start to be efficient in expelling dark matter from the 
central region with the final result of producing shallower 
density slopes. For 0.6 Gyr < t < 3.5 Gyr the value of a 
in the FAST simulation increases with time. On the con- 
trary, in the SLOW case, a only increases in the interval 
1.2 Gyr < t < 1.7 Gyr, and then remains approximately 
stable for t > 1.7 Gyr. This difference is caused by the 
different SMBH mass evolution observed in the two runs 
which makes AGN feedback operate in different modes, as 
already observed in Figure Q] Despite these differences, the 
final state of the system is very similar in both the FAST 
and SLOW cases, demonstrating that once AGN feedback 
is properly triggered and repeated accretion-ejection cycles 
are produced, the final configuration of the system is weakly 
dependent on the initial conditions. 



3.1 Experiments with smooth gas distributions 

To understand the importance of clumpy accretion, we also 
carried out simulations using a smooth halo without infalling 
gas clouds. These simulations showed that it is not possible 
to trigger heating-ejection cycles when low metallicity [Z ~ 
O.OIZq) is adopted: the cooling time is so long that the 
gas cannot cool efficiently. AGN feedback stays active in its 
weak radio mode, which provides enough energy to maintain 
equilibrium and the original gas distribution hardly evolves. 
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Gas Density (H/cmA3) -t = 2.55 Gyr 



Gas Density (H/cmA3) - 1 = 2.5S Syr 




Scaled M BH - SLOW 
M sa , [M ] - SLOW 




2.54 2.56 2.58 2.60 2.62 2.64 2.66 2.68 2.70 

t [Gyr] 



Figure 3. Projected density maps are shown for four snapshots over an accretion event followed by strong feedback. The SMBH is 
located at the centre of each map. This time sequence illustrates how AGN activity is able to excavate a hole in the gas density that is 
replenished on a time scale of ~ 10 s Gyr. The right panel shows the evolution of the gas mass within 5 kpc from the SMBH together 
with the gas accretion rate (rescaled). 



FAST 

SLOW 





Figure 4. Left panel: evolution of the central slope of the dark matter profile as a function of time when the dark matter density peak is 
chosen as a centre. Right panel: evolution of the central slope of the dark matter profile as a function of time when the SMBH is chosen 
as a centre. 
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r [kpc] t [Gyr] 

Figure 5. Left: dark matter density profiles at different ages and for the SMOOTH AGNOFF and SMOOTH AGNON simulations which 
do not include gas clumps. The peak of the dark matter distribution is chosen as a center. In the AGNOFF case we observe the effects of 
adiabatic contraction. In the AGNON case the effect of the heating-ejection mechanism is very weak compared to the case with clumps. 
Right: evolution of the dark (red line) and gaseous (green line) mass within 5 kpc from the cluster center in the SMOOTH AGNON 
simulation. 



The situation is different if a higher metallicity is taken for 
the gas. 



Figure [5] shows the results we obtain choosing Z = 
0.2Z Q in the SMOOTH AGNON and SMOOTH AGNOFF 
simulations. For this choice of gas metallicity, the cooling 
time is comparable to that of the FAST and SLOW simu- 
lations after the gas in the clumps mixed with the gaseous 
halo. In the SMOOTH AGNOFF case gas cools fast enough 
to produce contraction of the dark matter profile and steep- 
ening of the central cusp after 3.5 Gyr. When AGN feedback 
is turned on, there is a slight flattening of the initial central 
cusp after 3.5 Gyrs, and within ~ 3 kpc from the centre. The 
right panel of Figure [5] shows how the mass enclosed within 
5 kpc evolves in the SMOOTH AGNON run. Gas infall char- 
acterizes the first 1.5 Gyr of evolution until AGN activity 
becomes strong enough to prevent further infall. Since dark 
matter still dominates the central regions of the halo, no 
significant variation of the enclosed dark matter mass is ob- 
served as a response to gas motions. By t — 2 Gyr the black 
hole is massive enough to trigger some weak heating-ejection 
cycles. 



These results need to be carefully interpreted. The re- 
gion affected by the AGN in these "cooling flow" simu- 
lations, is close to our resolution limit. Our adopted halo 
properties are chosen to correspond to a high redshift clus- 
ter progenitor, rather than a massive cooling flow cluster. 
The metallicity for the above simulation, that shows a mild 
evolution, is therefore rather high. 



4 SUMMARY 

We analysed the results of N-body hydrodynamical simula- 
tions of a dark matter halo of mass M200 = 1.42 x 10 13 M© 
in the presence of gas cooling and AGN feedback. Our ini- 
tial conditions are chosen to represent the early evolution 
of the mass distribution in cluster sized halos. Progenitors 
of halos of mass ~ 10 14 M© identified at redshift z — 
have a typical mass of a few 10 13 M at redshift z > 2. A 
low metallicity gaseous halo is embedded in the dark matter 
distributions and we provide additional fuel for the central 
SMBH by letting ten high metallicity cold gas clumps fall 
into the central region. In all our runs, we consider radia- 
tive gas cooling. In two of our runs we also set a SMBH at 
the centre of the halo and we implement AGN feed back fol- 
lowing a m odification of th e Booth fc Schavel (|2009h already 
adopted by iTevssier et al.l (|201ll ). 

The aim of this work was to further investigate and 
confirm the results of the cosm ological hydrodynamical sim- 
ulations of (IMartizzi et al.l 120121 ) . in which the central dis- 
tribution of stars and dark matter within a cD galaxy were 
modified through violent gas motions induced by AGN feed- 
back. These results are similar to the those obtained by 
other authors using supernovae feedback in low mass dark 
matter halos ( Governato et al.l 120101: IPontzen fc Governatol 
120121 ; iMaccio' et al.ll201ll ; ITevssier et al.ll2012D . " 

There are limitations to our study. The resolution of 
our simulations is approximately 0.5 kpc which could be 
improved in the future. Furthermore, compared to the ob- 
servational data, our simulations over-predict the size of the 
regions in which the shallowing of dark matter and stellar 
mass profiles is measured. The simple AGN feedback model 
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we adopt is probably too efficient in unbinding baryons from 
the central regions of halos. Given these considerations, our 
results should be considered as a demonstration of the max- 
imum effect that AGN feedback can have on the mass dis- 
tribution in halos hosting galaxy clusters. 

Given these caveats, we can summarise our results as 
follows: 

• In the absence of AGN feedback, the gas clumps fall 
into the central region of the halo and mix with the rest 
of the gas. Given that the clumps provide high metallicity 
material, the result is a decrease of the cooling time in the 
central regions of the halo. Radiative cooling leads to the 
condensation of gas at the halo centre. The dark matter halo 
responds to the gas condensation via adiabatic contraction, 
increasing the central density. 

• In the presence of AGN feedback, after 1-2 Gyr we find 
a reduction in the dark matter density in the central regions 
of the halo. Roughly 60 % of the dark mass within 5 kpc 
from the center is moved to larger radii. The result is a mass 
distribution characterized by a dark matter density profile 
with a density slope in the region r < 5— 10 kpc significantly 
shallower than that of the initial NFW model. 

• The only process responsible for the dark mass deple- 
tion in the central regions of the halo is related to AGN feed- 
back and the resulting gas dynamics. After the gas clumps 
fall into the central region of the halo and mix with halo gas, 
we observe cycles of fast accretion onto the SMBH, followed 
by AGN feedback activity that blows gas away. The expelled 
gas cools and falls back into the central region of the halo. 
The whole cycle happens on a timescale comparable to the 
dynamical time. Repeated episodes generate a series of grav- 
itational potential fluctuations that i nfluence the total mass 
distri bution in an irreversible way (jPontzen fc Governatol 
12012ft . 

• If we use a smooth gas distribution without infalling gas 
clumps we are not able to trigger the mechanism we just de- 
scribed unless we assume high gas metallicity (shorter cool- 
ing time). The observed effect on the dark matter profile is 
then rather small. High central gas densities that arise from 
infall emerges as a key element to produce cusp-core trans- 
formations. At high redshift this may be triggered through 
mergers, gas stripping and fueling from cold streams. 



R ecent observa ti onal work (I Sand et al.l 
2«l0j; iNewman et"afl |2009| . l201ll ; iRichtler et al.l 
20121 ) 



Newman et al 



2004 



2011 



shows that cluster with shallow cen- 
tral dark matter density profiles may have been observed. 
Furthermore, very massive elliptical galaxies are observed 
to have surface b rightn e ss profile s with shallow central 



slope s (iKormendvl fl999l; lOuillen et al" | |2000|; lLaine et al 



2003| ; iTruiillo et al.l I200I; lLauer et al.| |2005l ; ICote et al 
iKoi 



2007 ; iKormei ldv et al.ll2009l ; lGrahamll201ll ) 



The two phenomena might be related, given that dark 
matter and stars can both be influenced by the same physical 
processes. Indeed, a strong prediction from our simulations 
is that the sizes of t he central cores in th e dark matter and 
the stars are similar (M artizzi et al.ll2012r ). Additional obser- 
vational evidence would be important to clarify this claim 
and to reveal possible correlations between AGN activity, 
stellar and dark matter cores and the global properties of 
clusters and the intra-cluster medium. 



ACKNOWLEDGMENTS 

The simulations were performed on the Schrodinger cluster 
at the University of Zurich. 



REFERENCES 

Booth C. M., Schaye J., 2009, Monthly Notices of the Royal 

Astronomical Society, 398, 53 
Bournaud F., Juneau S., Le Floc'h E., Mullaney J., Daddi 

E., Dekel A., Due P.-A., Elbaz D., Salmi F., Dickinson 

M., 2012, ApJ, 757, 81 
Cote P., Ferrarese L., Jordan A., Blakeslee J. P., Chen C- 

W., Infante L., Merritt D., Mei S., Peng E. W., Tonry 

J. L., West A. A., West M. J., 2007, ApJ, 671, 1456 
Dubois Y., Devriendt J., Slyz A., Teyssier R., 2012, MN- 

RAS, 420, 2662 
Duffy A. R., Schaye J., Kay S. T., Dalla Vecchia C, Battye 

R. A., Booth C. M., 2010, MNRAS, 405, 2161 
El-Zant A. A., Hoffman Y., Primack J., Combes F., Shlos- 

man I., 2004, ApJ, 607, L75 
Fromang S., Hennebelle P., Teyssier R., 2006, Astronomy 

and Astrophysics, 457, 371 
Gnedin O. Y., Zhao H., 2002, MNRAS, 333, 299 
Goerdt T., Moore B., Read J. I., Stadel J., 2010, ApJ, 725, 

1707 

Governato F., Brook C, Mayer L., Brooks A., Rhee G., 
Wadsley J., Jonsson P., Willman B., Stinson G., Quinn 
T., Madau P., 2010, Nature, 463, 203 

Graham A. W., 2011, ArXiv e-prints 

Komatsu E., Seljak U., 2001, MNRAS, 327, 1353 

Kormendy J., 1999, in D. R. Merritt, M. Valluri, & 
J. A. Sellwood ed., Galaxy Dynamics - A Rutgers Sympo- 
sium Vol. 182 of Astronomical Society of the Pacific Con- 
ference Series, The Central Structure of Elliptical Galaxies 
and the Stellar-Dynamical Search for Supermassive Black 
Holes, pp 124 — h 

Korrnendy J., Fisher D. B., Cornell M. E., Bender R., 2009, 
ApJS, 182, 216 

Krumholz M. R., McKee C. F., Klein R. I., 2004, The As- 
trophysical Journal, 611, 399 

Laine S., van der Marel R. P., Lauer T. R., Postman M., 
O'Dea C. P., Owen F. N., 2003, AJ, 125, 478 

Lauer T. R., Faber S. M., Gebhardt K., Richstone D., 
Tremaine S., Ajhar E. A., Aller M. C, Bender R., Dressier 
A., Filippenko A. V., Green R., Grillmair C. J., Ho L. C, 
Kormendy J., Magorrian J., Pinkney J., Siopis C, 2005, 
AJ, 129, 2138 

Maccio' A. V., Stinson G., Brook C. B., Wadsley J., Couch- 
man H. M. P., Shen S., Gibson B. K., Quinn T., 2011, 
ArXiv e-prints 

Martizzi D., Teyssier R., Moore B., 2012, MNRAS, 420, 
2859 

Martizzi D., Teyssier R., Moore B., Wentz T., 2012, MN- 
RAS, 422, 3081 

Navarro J. F., Eke V. R., Frenk C. S., 1996, MNRAS, 283, 
L72 

Newman A. B., Treu T., Ellis R. S., Sand D. J., 2011, ApJ, 
728, L39+ 

Newman A. B., Treu T., Ellis R. S., Sand D. J., Nipoti C, 
Richard J., Julio E., 2012, ArXiv e-prints 



Cusp-core transformations induced by AGN feedback 9 



Newman A. B., Treu T., Ellis R. S., Sand D. J., Richard 
J., Marshall P. J., Capak P., Miyazaki S., 2009, ApJ, 706, 
1078 

Pontzcn A., Governato F., 2012, MNRAS, 421, 3464 
Quillen A. C, Bower G. A., Stritzinger M., 2000, ApJS, 
128, 85 

Ragone-Figueroa C, Granato G. L., Abadi M. G., 2012, 

ArXiv e-prints 
Read J. I., Gilmore G., 2005, MNRAS, 356, 107 
Richtler T., Salinas R., Misgeld I., Hilker M., Hau G. K. T., 

Romanowsky A. J., Schuberth Y., Spolaor M., 2011, 

A&A, 531, A119+ 
Sand D. J., Treu T., Ellis R. S., Smith G. P., Kneib J.-P., 

2008, ApJ, 674, 711 
Sand D. J., Treu T., Smith G. P., Ellis R. S., 2004, ApJ, 

604, 88 

Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 
Teyssier R., 2002, Astronomy and Astrophysics, 385, 337 
Teyssier R., Fromang S., Dormy E., 2006, Journal of Com- 
putational Physics, 218, 44 
Teyssier R., Moore B., Martizzi D., Dubois Y., Mayer L., 

2011, MNRAS, 414, 195 
Teyssier R., Pontzen A., Dubois Y., Read J., 2012, ArXiv 
e-prints 

Toro E. F., Spruce M., Speares W., 1994, Shock Waves, 4, 
25 

Tremaine S., Gebhardt K., Bender R., Bower G., Dressier 
A., Faber S. M., Filippenko A. V., Green R., Grillmair 
C, Ho L. C, Kormendy J., Lauer T. R., Magorrian J., 
Pinkney J., Richstone D., 2002, ApJ, 574, 740 

Trujillo I., Erwin P., Asensio Ramos A., Graham A. W., 
2004, AJ, 127, 1917 



